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' Since the expected rates for neutralino-nucleus scattering are expected to be small, 

' one should exploit all the characteristic signatures of this reaction. Such are: (i) 

In the standard recoil measurements the modulation of the event rate due to the 
Earth's motion, (ii) In directional recoil experiments the correlation of the event 
f"*^ ■ rate with the sun's motion. One now has both modulation, which is much larger 

and depends not only on time, but on the direction of observation as well, and 
a large forward-backward asymmetry, (iii) In non recoil experiments gamma rays 
following the decay of excited states populated during the Nucleus-LSP collision. 
1 Branching ratios of about 6 percent are possible. 
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1. Introduction 

It is now established that dark matter constitutes about 30 % of the en- 



ergy matter in the universe. The evidence comes from the cosmological 
observations , which when combined lead to: 

fl b = 0.05, n C DM = 0.30, tt A = 0.65 

and the rotational curves 2 . It is only the direct detection of dark matter, 
which will unravel the nature of the constituents of dark matter. In fact one 
such experiment, the DAMA, has claimed the observation of such signals, 
which with better statistics has subsequently been interpreted as modula- 
tion signals 3 . These data, however, if they are due to the coherent process, 
are not consistent with other recent experiments, see e.g. EDELWEISS and 
CDMS 4 . 

Supersymmetry naturally provides candidates for the dark matter con- 
stituents. In the most favored scenario of supersymmetry the LSP can 
be simply described as a Majorana fermion (LSP or neutralino), a linear 
combination of the neutral components of the gauginos and higgsinos 5_8 . 
We are not going to address issues related to SUSY in this paper, since 
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thy have already been addressed by other contributors to these proceeding. 
Most models predict nucleon cross sections much smaller the the present 
experimental limit as < W~ 5 pb for the coherent process. As we shall see 
below the constraint on the spin cross-sections is less stringent. 

Since the neutralino is expected to be non relativistic with average ki- 
netic energy < T >~ 40KeV(m x /W0GeV), it can be directly detected 
mainly via the recoiling of a nucleus (A,Z) in elastic scattering. In some 
rare instances the low lying excited states may also be populated 9 . In this 
case one may observe the emitted 7 rays. 

In every case to extract from the data information about SUSY from the 
relevant nucleon cross section, one must know the relevant nuclear matrix 
elements 10 -n . The static spin matrix elements used in the present work 
can be found in the literature 9 . 

Anyway since the obtained rates are very low, one would like to be able 
to exploit the modulation of the event rates due to the earth's revolution 
around the sun 12 ' 13 . In order to accomplish this one adopts a folding 
procedure, i.e one has to assume some velocity distribution 12 > 13 > 14 > 15 for 
the LSP. One also would like to exploit other signatures expected to show 
up in directional experiments 16 . This is possible, since the sun is moving 
with relatively high velocity with respect to the center of the galaxy. 



2. Rates 

The differential non directional rate can be written as 

dRundir = t~~ da(u,v)\v\ (1) 
m x ArriN 

where da(u, v) was given above, p(0) — 0.3GeV/cm 3 is the LSP density in 
our vicinity, m is the detector mass and m x is the LSP mass 

The directional differential rate, in the direction e of the recoiling nu- 
cleus, is: 

p(0) m l s/u 
dUdir = \v\v.e <d(v.e) — da{u 7 v o{ -= — v.e) 



in 



Am N 2ir ' jj, r vV2 

where Q(x) is the Hcavisidc function and 

du 
2(/j, r bv) 

where u the energy transfer Q in dimensionless units given by 



d<r(u, v) == — — 2 \{t s F{u) 2 + S 5pin F n («)] (2) 



, Q = {m p Ab}- 2 = 40A- 4 / 3 MeV (3) 

Qo 
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with b is the nuclear (harmonic oscillator) size parameter. F(u) is the 
nuclear form factor and Fn(u) is the spin response function associated 
with the isovector channel. 

The scalar cross section is given by: 

I ; ■ Li i 

fjb r 



2„S 



1 + 4^ 



1 



(since the heavy quarks dominate the isovector contribution is negligible). 
CTjv x o is the LSP-nucleon scalar cross section. The spin Cross section is 
given by: 



= ( 



Mr ^2 spin 



fl r (p) 



! Cspin — 



P ' X 3(l + 4) 2 



S(«) 



(5) 



S(u) « 5(0) - [(^^o(O)) 2 + 2^^(0)^(0) + r>!(0)) 2 ] 



(6) 



f° 

'id 

, f l --<j\~/v 1 - J- 

/a /a 

/a' /A are the isoscalar and the isovector axial current couplings at the 
nucleon level obtained from the corresponding ones given by the SUSY 
models at the quark level, f%(q), f\(q), via renormalization coefficients g A , 
9a, Lc - fl = 9a/a (q)Ja = 9a/A(?)- These couplings and the associated 
nuclear matrix elements are normalized so that, for the proton at u = 0, 
yield £spm = 1- If the nuclear contribution comes predominantly from 
protons = Cl = ftp), S(u) w f2 2 and one can extract from the data the 
proton cross section. If the nuclear contribution comes predominantly from 
neutrons ((Qo = — ^1 = ^n) one can extract the neutron cross section. 
In many cases, however, one can have contributions from both protons 
and neutrons. The situation is then complicated, but it turns out that 
g A = 0.1 , g\ = 1.2 Thus the isoscalar amplitude is suppressed, i.e 
5(0) w fi 2 . Then the proton and the neutron spin cross sections are the 
same. 



3. Results 

To obtain the total rates one must fold with LSP velocity and integrate 
the above expressions over the energy transfer from Q m i n determined by 
the detector energy cutoff to Q ma x determined by the maximum LSP ve- 
locity (escape velocity, put in by hand in the Maxwellian distribution), i.e. 
v esc = 2.84 v with v the velocity of the sun around the center of the 
galaxy(229 Km/s). 
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3.1. Non directional rates 



Ignoring the motion of the Earth the total non directional rate is given by 



R = K 



C C oh(A M-4))o-p iX o 



+ c S pi n (A, fi r {A))a p P x0 



^sptn 



where K = £i°i jil JITpx and 

m x o m v V \ / 
~U r (A) 



c coh (A, u r (A)) = 



Hr(p) 



At C oh{A) , Cspin {A, \X T {A)") — 



U r (A) 



(7) 



(8) 

where t is the modification of the total rate due to the folding and nuclear 
structure effects. It depends on Q m i n , i.e. the energy transfer cutoff im- 
posed by the detector and a = [firbvoV^]^ 1 The parameters c co h(A, fi r (A)), 
c S pin (A, (j, r (A)), which give the relative merit for the coherent and the spin 
contributions in the case of a nuclear target compared to those of the pro- 
ton, are tabulated in table 1 for energy cutoff Q m i n = 0, 10 keV. Via Eq. 
(7) we can extract the nucleon cross section from the data. 

Using ill = I- 22 an d ^i = 2 -8 for m I and 19 F respectively the extracted 
nucleon cross sections satisfy: 



p,x u 



oh (A,n r (A)) 

nn (A, Mr {A)) 



3 

Of 



-2 ^« xlO 4 (A = 127) , w xlO 2 (A = 19) 



(9) 

It is for this reason that the limit on the spin proton cross section extracted 
from both targets is much poorer. For heavy LSP, > 100 GeV, due to the 
nuclear form factor, t sp i n (l27) < i sp j n (19). This disadvantage cannot be 
overcome by the larger reduced mass (see Fig. 1). It even becomes worse, 
if the effect of the spin ME is included. For the coherent process, however, 
the light nucleus is no match (see Table 1 and Fig. 2 ). 



3.2. Modulated Rates. 

If the effects of the motion of the Earth around the sun are included, the 
total non directional rate is given by 

R = K [c coh {A, a r {A))a^ x0 {\ + h(a, Q mi n)cosa)] (10) 

and an analogous one for the spin contribution, h is the modulation am- 
plitude and a is the phase of the Earth, which is zero around June 2nd. 
The modulation amplitude would be an excellent signal in discriminating 
against background, but unfortunately it is very small, less than two per 
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Table 1. The factors cl9 = c coh (19, Mr(19)), sl9 = c spln (19, /i r (19)) and 
cl27 = c coh (127, ^t r (127)), sl27 = c spin (127, £t r (127)) for two values of Q = Q min . 



Q 
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cl9 


2080 


2943 


3589 


4083 


4471 


5037 


5428 


6360 





sl9 


5.7 


8.0 
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10.9 


11.9 


13.4 


14.4 


16.7 





cl27 


37294 


63142 
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114295 


131580 


142290 
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sl27 


2.2 
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Figure 1. The coefficient c ap i„ (A, fi r (A)) as a function of the LSP mass. The dashed 
curves correspond to the 19 F system ( the short for Q m i n = 0, and the long for Q m i„ = 
10 keV). The solid curves correspnd to the 127 I (the thin for Q m i n = and the thick for 

Qmin = 10 keV. 



cent. Furthermore for intermediate and heavy nuclei, it can even change 
sign for sufficiently heavy LSP (see Fig. 3). So in our opinion a better sig- 
nature is provided by directional experiments, which measure the direction 
of the recoiling nucleus. 
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Figure 3. The modulation amplitude h as a function of the LSP mass in the case of 
127 I for Qmin = on the left and Q m i n = 10 kcV on the right. 



3.3. Directional Rates. 

Since the sun is moving around the galaxy in a directional experiment, i.e. 
one in which the direction of the recoiling nucleus is observed, one expects 
a strong correlation of the event rate with the motion of the sun. In fact 
the directional rate can be written as: 

Rdir = [ccoh{A, fj, r (A))(Tp iX o(l + h m cos(a - a m ir))] (11) 

and an analogous one for the spin contribution. The modulation now is h m , 
with a shift a m ir in the phase of the Earth a, depending on the direction of 
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Table 2. The parameters t, h, re, h m and a m for Q m i n = 0. 
The results shown are for the light systems. +x is radially out 
of the galaxy , +z is in the the sun's motion and +y vertical 
to the plane of the galaxy so that ((x,y,x) is right-handed. 
C-m = 0, 1/2, 1, 3/2 implies that the maximum occurs on June, 
September, December and March 2nd respectively. 



type 


t h dir 


K 








+z 


0.0068 


0.227 


i 


dir 


+ (")x 


0.080 


0.272 


3/2(1/2) 




+(-)y 


0.080 


0.210 


0(1) 




-z 


0.395 


0.060 





all 
all 


1.00 

0.02 









observation. k/(2tt) is the reduction factor of the unmodulated directional 
rate relative to the non-directional one. The parameters k , h m , a m 
strongly depend on the direction of observation. We prefer to use the 




0.5 1 1.5 2 2.5 3 0.5 1 1.5 2 2.5 3 



Figure 4. The expected modulation amplitude h m for A = 127 in a direction outward 
from the galaxy on the left and perpendicular to the galaxy on the right as a function 
of the polar angle measured from the sun's velocity. For angles than n/2 it is irrelevant 
since the event rates are tiny. 

parameters k and h m , since, being ratios, are expected to be less dependent 
on the parameters of the theory. In the case of A = 127 we exhibit the the 
angular dependence of h m for an LSP mass of m x = lOOGeV" in Fig. 4. 
We also exhibit the parameters t, h, k, h m and a m for the target A = 19 
in Table 2 (for the other light systems the results are almost identical). 
The asymmetry in the direction of the sun's motion 17 is quite large, w 0.97, 
while in the perpendicular plane the asymmetry equals the modulation. 
For a heavier nucleus the situation is a bit complicated. Now the parameters 
k and h m depend on the LSP mass as well, (see Figs 5 and 6). The 
asymmetry and the shift in the phase of the Earth are similar to those of 
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the A = 19 system. 



0.065 



0.055 




Figure 5. The parameter k as a function of the LSP mass in the case of the A = 127 
system, for Q m i n = expected in a plane perpendicular to the sun's velocity on the left 
and opposite to the sun's velocity on the right. 



0.19 
0.18 




Figure 6. The modulation amplitude h m in a plane perpendicular to the sun's velocity 
on the left and opposite to the suns velocity on the right. Otherwise the notation is the 
same as in Fig 5. 



3.4. Transitions to excited states 

Incorporating the relevant kinematics and integrating the differential event 
rate dR/du from u m in to u max we obtain the total rate as follows: 

ru max ,n 2 ru max in 

*■— /_ J t (i -^ )du - R >-=L ,12) 

where u exc — j^^q and E x is the excitation energy of the final nucleus, 
u m ax = (y/a) 2 - (E x /Qo) , V = v/upsilono and u min = Q mm /Qo, Qmin 
(imposed by the detector energy cutoff) and u max = {y eS c/o) 2 is imposed 
by the escape velocity (y esc = 2.84). 

For our purposes it is adequate to estimate the ratio of the rate to the 
excited state divided by that to the ground state (branching ratio) as a 
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function of the LSP mass. This can be cast in the form: 

'S'exc(O) ^ exc{^excj ^umax) [1 ~t" h exc {u eX c> u max) COS Q;] 



BRR 



S gs (0) 



* j S («mm)[l + h(u m in) COS tt] 



(13) 



in an obvious notation ). S gs (0) and S exc (0) are the static spin matrix 
elements As we have seen their ratio is essentially independent of super- 
symmetry, if the isoscalar contribution is neglected. For 127 I it was found 
9 to be be about 2. The functions are given as follows : 



^ gs (^min) 



^ exc{ u exc: ^max) — 



S gs (u) 
ma " S exc {u) 



F?°(u) ty(av^) - MVesc)] du (14) 



s (0) ^iT c W(l-^f )WMl+u exc /u))-^{y esc )] du 

(15) 



The functions ip arise from the convolution with LSP velocity distribution. 
The obtained results are shown in Fig. 7. 




m LS p -> (GeV) 

Figure 7. The ratio of the rate to the excited state divided by that of the ground state 
as a function of the LSP mass (in GeV) for 127 1. We assumed that the static spin matrix 
element of the transition from the ground to the excited state is a factor of 1.9 larger 
than that involving the ground state, but the functions Fn(ti) are the same. On the left 
we show the results for Q m ;„ = and on the right for Q m i n = 10 KeV. 



4. Conclusions 

Since the expected event rates for direct neutralino detection are very 
low 5 ' 8 , in the present work we looked for characteristic experimental signa- 
tures for background reduction, such as: 

Standard recoil experiments. Here the relevant parameters are t and h. 
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For light targets they are essentially independent of the LSP mass 17 , essen- 
tially the same for both the coherent and the spin modes. The modulation 
is small, h w 0.2%, but it may increase as Q m i n increases. Unfortunately, 
for heavy targets even the sign of h is uncertain for Q m i n = 0. The sit- 
uation improves as Q m in increases, but at the expense of the number of 
counts. 

Directional experiments 16 . Here we find a correlation of the rates with 
the velocity of the sun as well as that of the Earth. One encounters re- 
duction factors k/2tt, which depend on the angle of observation. The most 
favorable factor is small, ~ 1/Att and occurs when the nucleus is recoiling 
opposite to the direction of motion of the sun. As a bonus one gets modu- 
lation, which is three times larger, h m =w 0.06. In a plane perpendicular 
to the sun's direction of motion the reduction factor is close to 1/127T, but 
now the modulation can be quite high, h m w 0.3, and exhibits very inter- 
esting time dependent pattern (see Table 2. Further interesting features 
may appear in the case of non standard velocity distributions 15 . 
Transitions to Excited states. We find that branching ratios for tran- 
sitions to the first excited state of 127 I is relatively high, about 10%. The 
modulation in this case is much larger h exc w 0.6. We hope that such a 
branching ratio will encourage future experiments to search for characteris- 
tic 7 rays rather than recoils. Acknowledgments: This work was supported 
in part by the European Union under the contracts RTN No HPRN-CT- 
2000-00148 and MRTN-CT-2004-503369. Part of this work was performed 
in LANL. The author is indebted to Dr Dan Strottman for his support and 
hospitality. 
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